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Abstract 

In its simplest form the curvaton paradigm requires the Hubble parameter during 
inflation to be bigger than 10 8 GeV, but this bound may be evaded in non-standard 
settings. In the heavy curvaton scenario the curvaton mass increases significantly 
after the end of inflation. We reanalyze the bound in this set up, taking into account 
the upper bound on the curvaton mass from direct decay. We obtain if* > 10 s GeV 
if the mass increase occurs at the end of inflation, and > 10 -14 GeV if it occurs 
just before nucleosynthesis. We then discuss two implementations of the heavy 
curvaton. Parameters are constrained in these explicit models, and as a result even 
obtaining TeV scale inflation is hard to obtain. 



1 Introduction 



In the curvaton scenario not the inflaton but some other light field — the curvaton 
field — is responsible for the observed density perturbations [I]. In the post infla- 
tionary epoch the curvaton oscillates around the origin of its potential in a radiation 
dominated background. It is during this period that the curvaton contribution to the 
total energy density grows and the isocurvature perturbations of the curvaton field 
are converted into adiabatic ones. The resulting perturbation spectrum is Gaussian 
if the curvaton accounts for at least one percent of the total energy density at the 
time of decay. 

The great merit of the curvaton scenario is that it liberates models of inflation j2] . 
Both the scale of inflation, and the flatness of the inflaton potential are decoupled 
from the observed magnitude and scale dependence of the cosmic microwave back- 
ground (CMB) spectrum. This means in particular that inflation can take place at 
a much lower scale than the customary if* ~ 10 10 — 10 14 GeV. 

Many models of inflation have been constructed, in which the inflaton potential 
is generated by the same supersymmetry (SUSY) breaking mechanism that oper- 
ates in the vacuum, and the scale of inflation is naturally ii* ~ 10 3 GeV jSl IH E]- 
However, as was shown in Ref. jH] the curvaton scenario in its simplest form cannot 
accommodate such low inflationary scales either. The reason is that the lower the 
scale of inflation, the lower is the initial curvaton energy density and the shorter 
the period of oscillations. Only for large enough Hubble constants during inflation, 
ii* > 10 7 GeV, can the curvaton come to dominate the energy density (and decay) 
before the epoch of nucleosynthesis. 

To rescue low scale models of inflation it is necessary to go beyond the simplest 
curvaton scenario. 1 One way to reconcile low scale inflation with the curvaton 
scenario is to consider curvaton oscillations in a background dominated by a stiff fluid 
instead of radiation jHj . In Refs. (HI E] the "heavy curvaton" scenario was proposed, 
in which the curvaton mass increases significantly after the end of inflation. In this 
paper we reanalyze the bound on the Hubble constant during inflation within this 
scenario. As the curvaton mass increases so does its decay rate. This leads to an 
upper bound on the curvaton mass at the time of curvaton decay; increasing the 
mass beyond this bound will lead to a rapid decay of the curvaton condensate. This 
constraint was not taken into account in previous works. We obtain the bound 
on the inflationary scale in the heavy curvaton scenario if* > 10 8 GeV if the mass 
increase occurs right after the end of inflation, and if* > 10~ 14 GeV if it occurs just 
before big bang nucleosynthesis (BBN). 

In the next section we start with a review of the bound on the Hubble scale 
during inflation in the simplest curvaton scenario. How this bound can be evaded 
is discussed in turn in section |3J In particular, the heavy curvaton scenario is 
detailed. The rest of the paper is devoted to two specific implementations of the 

1 See Ref. (7 4 for a discussion of low scale inflation in the context of the inhomogeneous reheating 
scenario. 
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heavy curvaton. In section|3]we discuss a model in which the curvaton mass increases 
due to a direct coupling to a "Higgs" field. In section we analyze a scenario in 
which the curvaton is an axion. Its mass receives contributions from two different 
non-renormalizable operators; one operator dominates during inflation, whereas the 
other dominates afterwards. Although theoretically an inflationary scale as low as 
~ 10~ 14 GeV is possible, we find that in these explicit models even low scale 
inflation with ~ TeV is hard to obtain. The reason is that once an explicit 
model is specified, there are relations between the various parameters, which can no 
longer be varied independently. 



2 The bound on the inflationary scale 

We define the curvaton density parameter Q a = p a j p tot , which gives the contribution 
of the curvaton to the total energy density. The curvaton density parameter changes 
with time; at the onset of curvaton oscillations 

a = mla * a) 

"osc m M 2 > V 1 ) 
- n osc JW pl 

with m a the curvaton mass. The subscripts 'osc' and '*' denote the corresponding 
quantity at the onset of curvaton oscillations, and during inflation at the time ob- 
servable scales leave the horizon, respectively. We have dropped the subscript a from 
the curvaton density parameter to avoid notational cluttering. If the curvaton and 
the background energy density red shift at different rates, the curvaton contribution 
to the total energy density changes, according to 

*t = -a(fO— , (2) 

with 

(tt\ 2 (a; tot — a; CT ) 

MH) = l+u,„ ■ (3) 

Here u tot = (EjCjjPi) / p to t with pi = ujipi the equation of state for each fluid. The 
quantity a(H) changes if either the equation of state parameter of the background or 
that of the curvaton changes. In particular, when the curvaton energy density comes 
to dominate the total energy density, u tot — > u a , the back ground and curvaton red 
shift at the same rate and a(H) — > 0. The ratio Q a remains less than unity as it 
should be. 

Consider the simple case that a remains constant after the onset of curvaton 
oscillations until the curvaton energy density approaches the total energy density. 
The generalization to a changing a, which occurs for example if the inflaton field 
decays after the onset of oscillations, is straightforward; we will comment on it later. 
Denote by a the constant value of a(H) when p a ^ p tot . Since a(H) vanishes in the 
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limit that the curvaton comes to dominate the energy density we can write down 
the inequality: 2 

mjal (H osc \ a 

dec ~CM p 2 ArJ ' [) 

with To- the curvaton decay rate, and Qdec the curvaton density parameter at the 
time of curvaton decay. The observed spectrum of density perturbations is obtained 

for mmm 

a* « Q^qA^H^, (5) 

with A = 1.57T x 10 -4 . The damping factor q is defined through (Sa/a)^ — 
q{8a/a) if . For a quadratic or flat potential q = 1, but it is different for other mono- 
mials. Combining Eqs. (jU |3J) leads to the bound on the Hubble constant during 
inflation: 

H^ HoscM ff ( Ta Y /2 ■ (6) 
m CT gfi^ \H osc J 

The minimum value for the decay rate is IV > max[r grav , Tbbn], with r grav ~ m^/M^ 
the decay rate for gravitational decay, and Tbbn ~ 4x 10~ 25 GeV corresponds to decay 
just before BBN. 

The bound on if* is minimized in the limit that the curvaton decays just when 
it comes to dominate the energy density. The period of curvaton oscillations has to 
be sufficiently long for domination to be possible. Since for a lower scale of inflation 
the initial curvaton vacuum expectation value (VEV) is smaller, see Eq. (JSJ), this 
period becomes increasingly long in this limit. It might seem that the period of 
curvaton oscillations can be shortened, and therefore the bound be lowered, if the 
curvaton decays when it is still sub-dominant in energy and Qdcc < 1- However, this 
advantage is more than off set by the fact that to obtain the right size of density 
perturbations the initial curvaton VEV is lower, and therefore Q osc is lower. As a 
result the bound is actually stronger by a factor fi^/ 2 . 

A Gaussian perturbation spectrum, in agreement with observations, requires 
fldec > 10~ 2 , that is, the curvaton has to account for at least one percent of the 
total energy density at the time of decay [TUj . 



2.1 Standard scenario 

In the simplest models the curvaton field remains frozen until H ~ m a , at which 
point the field starts oscillating in the potential well. The curvaton potential is 
quadratic and q ~ 1. After inflaton decay the background is radiation dominated 

2 If r CT is much smaller than -ffdom, the Hubble constant at which the curvaton comes to dominate 
the energy density, this equality is satisfied by many orders of magnitude. As a consequence the 
Hubble bound in Eq. |JBJ is off (to weak) by an amount {T (J /H l \ - m ) a / 2 . We will come back to this 
issue when discussing the heavy curvaton. 
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and a = 1/2. We will refer to these parameter values as the standard scenario. The 
bound in Eq. JBJ then becomes 

H* > AM pl \ 1 ' 2 . (7) 

Here T a ~ X 2 m a and A > max [Abbn, A grav ], with 

Abbn ~ V H bbn/m a , (8) 
Agrav ~ m a /M ph (9) 

the coupling corresponding to decay just before BBN, and a gravitational strength 
coupling respectively. Plugging in the minimum value for the coupling gives 



H* > max 



;i(rc; (A (^y /5 , 5 xlO n GeV(^ 



(10) 



The first bound in Eq. (fTUj) is stronger for Abbn > A grav or m a < m cr (ifbbn) ~ 
10 4 GeV, whereas in the opposite limit the second bound applies. Here m CT {T a ) is 
the mass for which the decay rate becomes of gravitational strength: 

m cr (r CT ) = (rvMp 2 ,) 173 . (ii) 

The bound is minimized by minimizing the coupling. For A = Abbn this means 
increasing the mass to its largest possible value; for A = A grav it means decreasing 
the curvaton mass as much as possible. Hence, the bound is weakest in the limit 
m a — > m CT (H hhn ) for which r grav — > T hhn . The absolute lower bound then is 

#*>8xl0 7 GeV. (12) 

This result is clearly incompatible with low scale inflation with ~ TeV. 

The lower bound is obtained in the limit Q^ec 1 ; when Eq. (@J) is saturated. 
It is assumed that inflaton decay happens before the onset of curvaton oscillations, 
Tj > H osc , with / denoting the inflaton field. Otherwise the lower bound is increased 
by a factor (m^/IV) 1 / 4 . 



3 Evading the bound 

To evade, or at least, to weaken the above bound on the Hubble constant during 
inflation, one can try one or several of the following possibilities: 

• q > 1 

• a > 1/2 

• H osc /m a < 1 
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A damping factor q > 1 is possible with a potential that grows faster than 
quadratic at large VEV ^2]- Nevertheless, the bound will be stronger as the com- 
bination cr* oc Q a q always decreases faster than for a quadratic potential. Another 
possibility is that the potential increases slower than quadratic. However, small sub- 
horizon scale fluctuations will grow at a much faster rate than the super-horizon scale 
fluctuations observed by CMB experiments. The condensate will fragment or decay 
into Q-balls when H ~ m a ^3]. Q-balls behave as non-relativistic matter, and the 
curvaton bound is not much different from the standard case. 

It is possible to have a > 1/2 if the background is dominated by a fluid that 
red shifts faster than radiation and u tot > 1/3. If the inflaton energy density is 
dominated by kinetic energy the inflaton behaves as a stiff fluid and u tot = 1, giving 
an explicit realization. 

If the effective curvaton mass increases after the end of inflation it is possible to 
have H osc /m a < 1. Such a mass change could be triggered by a phase transition - 
a Higgs mechanism comes to mind. Thermal masses are no good for this purpose, 
as early thermal evaporation should be avoided. 



3.1 Stiff inflaton fluid 

The bound in Eq. (jUJ) is weakened by increasing the exponent a > 1/2. The max- 
imum value, a = 1, is obtained for a stiff fluid with uo = 1. This is exactly what 
happens if inflation is followed by a period of kination, which occurs for example in 
quintessential inflation G2] and disformal inflation [5] . The bound for a back- 
ground dominated by a stiff fluid is 

> M pl AX. (13) 
Plugging in the minimal coupling gives 



H* > max 



H \ 1/3 

10 2 GeV ( — , 5 x KT 4 m ff 



(14) 



For m a < m cr (if bbn ) ~ 10 4 GeV (see Eq. flTTj) ) the first bound is stronger, whereas 
in the opposite limit the second expression prevails. A scale invariant spectrum is 
obtained for > lOm^, which is actually a stronger bound for gravitational decay. 
The absolute lower bound on the scale of inflation then is 

> 10 2 GeV. (15) 

This is compatible with low scale inflation, though not very comfortably. 

In Ref. jS] the bound on the Hubble scale for gravitational decay — the second 
bound on the right hand side of Eq. (|3.1J) — was derived. However, Ref. jH| then 
proceeds by taking m — > m cr (ifbbn) ~ TeV to obtain the lower bound > GeV, 
forgetting that m ^> if* is incompattible with a scale invariant perturbation spec- 
trum. 
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It is assumed that inflaton decay happens after the curvaton comes to dominate 
the energy density. Otherwise, if the inflaton decays into radiation, a — > 1/2 as in 
the standard scenario and the bound gets stronger by a factor (I^/TV) 1 / 4 . For TeV 
scale inflaton masses and gravitational strength couplings inflaton decay happens at 
Tj ~ 10~ 4 ift>bn- This poses no problem as long as the inflaton energy density is small 
enough at the time of decay, so that the accompanying entropy production is small; 
in this scenario the universe is reheated by curvaton decay. Low scale inflation with 

~ TeV requires T a /H hhn < 10 2 . This constitutes tuning. 

3.2 Heavy curvaton 

In the heavy curvaton scenario the curvaton mass increases significantly during a 
phase transition, at H = H pt , from (mj), — > m a ^> (iijj), with [m a \ the curvaton 
mass during inflation. If the curvaton mass after the phase transition is smaller 
than the Hubble constant m CT < -ffpx, oscillations set in when H osc ~ m a and the 
bound of the standard scenario applies, with the sole difference that now the mass 
can be larger than the Hubble scale during inflation (but {m a \ < to assure a 
scale invariant spectrum). 

The more interesting case occurs in the opposite limit, m a > H pt , when oscilla- 
tions set in at the time of the phase transition and the ratio H osc /m a is less than 
unity. This may weaken the bound on the inflationary scale Eq. © with respect to 
the standard scenario. This possibility was discussed in Refs. jHIEl- However, these 
papers did not take into account the accompanying increase in the decay rate, and 
the constraints this puts on the scenario. 

3.2.1 Mass increase at the end of inflation: H pt ~ 

We consider first the case that the curvaton mass increases right after the end of 
inflation. The advantage of such a scenario is that the phase transition triggering 
the mass increase can be be part of the inflaton sector, and hence no new dynamics 
need to be introduced. This would be realized for example if the curvaton field is 
coupled to the waterfall field in hybrid inflation. 

As the curvaton mass increases so does the curvaton decay rate. As follows from 
Eq. ©, a larger mass can lower the bound on H*, whereas a larger decay rate is 
counter productive. The upper bound on the mass is 



m a < mm[m en , m& 



(16) 



with 



171, 



en 



AH pt M pl 



(17) 



Ice 
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The first bound, m a < m en , is that the energy density in the curvaton should be less 
than the total energy in the universe at the time of the phase transition, Q pt < 1 
with 



J 'pt — 

Ptot 



H=H vt V H pt M p i 



(19) 



If the decay rate increases above r CT > H pt curvaton decay follows immediately; this 
sets the upper bound on the curvaton mass at the time of decay m a < m dcc . The 
first expression on the right hand side of Eq. (|18|) corresponds to a coupling constant 
A > A grav , whereas the second expression is for gravitational strength couplings. 

The bound on the inflationary scale is minimized for the largest possible curvaton 
mass. The mass md ec is maximized in the limit A — > A grav . Direct decay with Q^ec ~ 
Q pt ~ 1 is only possible if the curvaton dominates the energy density and m dec ~ 
m en . With H osc ~ T a ~ H* and standard scenario values for other parameters the 
bound on the Hubble scale during inflation is: 

> A 3 M pl = 2 x 10 8 GeV. (20) 

This can be found using Eq. (jHJ) together with the upper bound on m a given by 
Eq. (dnj). Alternatively, one can solve md ec ~ m en with A — > A grav . For smaller 
values of the Hubble constant m& cc < m en and decay will happen before the curvaton 
dominates the energy density. 

The above expression is valid for fl^ ec ~ 1. For smaller values of fl^ec < 1 

— 1/2 

the mass constraint from energy conservation is relieved m en cx f2 dec ' , while m,x ec 
remains unchanged. As a result, md ec < men and decay occurs when the curvaton 
density parameter is too small Q a < Qdec, unless 

A 3 M pl , . 

* > ~^f- (21) 

^dec 

The bound on the inflationary scale is minimized in the limit Q^ec — > 1, and is given 
by Eq. (|2TH). Our results differ from Lyth [Hj, who finds > 10 5 GeV. The reason is 
that although Lyth considers energy conservation (he takes fl^ec —> 10~ 2 to increase 
?n en as much as possible), he does not take into account the constrains from direct 
decay, i.e., the constraint m < m^ec- 

Low scale inflation is excluded. What is more, this scenario offers no improve- 
ment over the standard scenario, while complicating the curvaton potential consid- 
erably. 

3.2.2 Mass increase after inflation: H pt H* 

From Eq. © it can be seen that the bound on the Hubble constant can be lowered 
if 

H osc = max[if pt , (m^)*] (22) 
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is lowered, i.e., if H pt <C H*. This was noticed in Ref. j^]. Delaying the onset of 
oscillations means lowering m^ ec , see Eq. (|18|). which is counter productive. However, 
the effect on H osc wins, and it is advantageous to lower H osc as much as possible. 
There is no lower bound on the scale of the phase transition except that it should 
take place before curvaton decay, which in turn should occur before BBN: H pt > 

The bound on the Hubble scale during inflation is minimized in the limit H osc = 
H pt —>■ ifbbn; and direct decay with f2d ec ~ ^ P t ~ 1- After minimizing the onset of 
oscillation, the mass md ec should be maximized, which means A — > A grav . (In this 
limit Abbn = A grav ). Hence, we find just as in the previous subsection that the bound 
is minimized for A ~ A grav and m en ~ m^ ec . This last requirement is equivalent to 
demanding the curvaton to decay right after the phase transition while dominating 
the energy density. In addition, and in contrast with the scenario in which the 
phase transition takes place promptly at the end of inflation, the bound is further 
minimized by taking H osc —>■ H^ n . The bound on the Hubble constant then is 



i > max[A 2 AM p i, AH^M^ 3 } > 10~ 14 GeV. (23) 



The equation is saturated for direct decay. The first bound is strongest if the 
coupling A > Ag rav or m a < m CT (T a ), and vice versa for the second bound. 

Low scale inflation with ii* ~ TeV is clearly possible. It requires H osc < GeV, 
and thus (m^)*, H pt < GeV. If {m a ) if > H pt , and oscillations set in before the 
phase transition, the bound is raised by a factor ((m^)*/ H pt ). 

If the energy density in the curvaton field after the phase transition is less than 
the total energy density, Q pt <C 1, direct decay will lead to an unacceptable level 
of non-Gaussianity. A period of curvaton oscillations is required, long enough for 
the curvaton to dominate at decay. By energy conservation, the curvaton mass is 
vn a = f2p{ 2 m en , with m en given in Eq. (|17|) . In a radiation dominated universe with 
a = 1/2 the curvaton comes to dominate the energy density at H^om = ^pt-^pt- At 
that moment a — > 0, and there is no further gain in lowering IV < H^ om . The lower 
bound on the Hubble scale during inflation, see Eq. ©, is oc H^^/m^ and thus 
the fl pt dependence cancels out. The bound on the Hubble constant for Q pt < 1 is 
the same as for direct decay with Q pt = 1, as long as domination takes place before 
decay, that is T a < ifdom = ^pt-^pt- This latter requirement is actually stronger and 
gives 

' AM p iX 2 AHll'M 1 /' 



if* > max 



3/2 ' l/6 
1 J pt 1 'pt 



(24) 



for a coupling A > Abbn, A grav and for A = A grav respectively. 

Having Q pt < 1 increases the bound for a gravitational decay rate. The reason is 
that r grav oc Q p { 2 ml n , and unless the Hubble constant during inflation is higher than 
that for direct decay (the value when Eq. (j2Hj) is saturated) decay happens before 
domination. Likewise, for a fixed coupling constant X > A grav , the bound is minimized 
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for direct decay with Cl pt ~ 1. If fi pt < 1, then since IV oc f2 p { m en , also here decay 
will happen before domination unless H* is raised above the saturation value in 
Eq. (J2HJ). Another way to interpret the bound Eq. (J2IJ), however, is that for a fixed 



inflationary scale H*, the curvaton scenario can equally work for f2 pt ~ 1 and A < 
H*/(AM p i), as for fi pt < 1 and a smaller coupling A 2 < fipfiJ,/(AM p i). In both 
cases the curvaton dominates before decay. The requirement max[if bbn , r grav ] < 
r CT < Hdom reads in terms of the coupling 



max 



-£f bbrl -£/* 



nH 2 AH pt M pl 



< A 2 < 



q 3//2 /7 



M pl A 



(25) 



The minimum bound Eq. is obtained in the limit -^bbn ~ Tgrav ~ H t 
requires tuning Vt pt = 1 and A = H*/(M pl A) = (H pt / 'Mpi) 1 / 3 . 

Eq. (|25|) gives a lower bound on Q pt for the curvaton scenario to work 



dom 



which 



f2 pt > max 



bbn 



P t 



1/2 



Mpiffp 2 t A 3 
# 3 



For low scale inflation ~ 10 3 GeV, together with Eq. 



(26) 

this gives fi pt > 10~ 12 . 



3.3 Additional ingredients 

Above we have shown that low scale inflation is compatible with the curvaton sce- 
nario in certain circumstances. But there are obstacles every working model has to 
face, such as the existence of one or more small masses and late time baryogenesis. 
Small couplings are needed to avoid thermal evaporation. 

A small curvaton mass is needed, (mo-)* < O.lif* for the stiff inflaton fluid, and 
(trier)* < 1 GeV in the heavy curvaton scenario. This is hard to obtain in the presence 
of supersymmetry/supergravity, as all scalar fields obtain a soft mass m ~ m 3 / 2 from 
low energy SUSY breaking, and m ~ from SUSY breaking by the finite energy 
density during inflation. Possible ways out are tuning (slightly in the case of a 
stiff inflaton fluid but considerably for the heavy curvaton); special, non-minimal 
Kahler potentials which suppress soft masses JH]; global symmetries (the curvaton 
as pseudo Nambu Goldstone boson) [T7j; giving up on supersymmetry. Soft masses 
can be small in gauge mediated SUSY breaking schemes, e.g. SUSY breaking at a 
scale F > (10 5 GeV) 2 yields a gravitino mass rrt.3/2 > 10~ 9 GeV. However, one can 
no longer associate the inflaton sector with the SUSY breaking sector, or employ 
moduli fields to obtain inflation, as this would also imply a very small inflationary 
scale H m ~ F/M ph 

Baryons and dark matter cannot be created — the epoch of creation being the 
epoch after which their comoving number density remains constant — when the 
curvaton contributes still little to the total energy density C 1, as this would 
give rise to an unacceptable level of isocurvature perturbations [TO! ITH]. Creation 
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by the curvaton field is only consistent with an adiabatic perturbation spectrum, if 
the curvaton is close to dominating the energy density at decay: Qd ec > 0.6 (0.9) for 
baryons (cold dark matter) 3 . There are no constraints if CDM or baryon production 
takes place after curvaton decay. This is a serious obstacle, especially for small 
curvaton decay rates and late time phase transitions T a , H pt — > H^ n , as there are 
no known production mechanisms that work at such a late time. 

The stiff inflaton fluid is not affected by thermal effects, as the reheating temper- 
ature is low and all couplings involved are small. This is not necessarily the case for 
the heavy curvaton. Consider a curvaton coupling A to a standard model (SM) field. 
The SM field is in thermal equilibrium with the radiation bath if T > Xa. Evapora- 
tion of the curvaton condensate is avoided if the scattering rate is sufficiently small: 
To- ~ \a g T < H (with a g a SM structure constant) [TH], or 

H > X 2 a 2 g M pl . (27) 

Note that if this equation is satisfied comfortably, even the case of direct decay is 
compromised, as the thermal bath from curvaton decay will give rise to evaporation 
before fl a ~ 1. 

4 The Higgsed curvaton 

A straightforward way to implement the heavy curvaton is by introducing a coupling 
h 2 <fr 2 a 2 in the potential. The field has zero VEV during inflation, but acquires a 
large VEV during some phase transition at H = H pt < if*. The effective mass of 
the curvaton increases suddenly at this phase transition. 

4.1 The potential 

Consider the following superpotential 

W = KS((j) 2 - v 2 ) + ha(<j) 2 - f). (28) 

We take v = f so that the energy density vanishes in vacuum; this will not influence 
the results in an essential way. The F-term potential is 

V = {\k\ 2 + \h\ 2 )\(f) 2 - v 2 \ 2 + A\(j)\ 2 \KS + ha\ 2 . (29) 

The potential is of the kind used in hybrid inflation. In analogy we will call S the 
flaton (flat direction) field, and <fi the waterfall field. The curvaton is denoted by a 
as before. 

3 Barring the possibility that baryon production before decay leads to an isocurvature perturba- 
tion which is exactly canceled by an isocurvature perturbation from CDM production after decay, 
or vice versa. 
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S and a are massless in the supersymmetric limit. They have small masses 
induced by Planck suppressed operators and/or soft masses, which we assume are 
much smaller than the scale of inflation. During inflation S and a are flat and 
acquire a large VEV. The field is heavy due to its coupling with S and remains 
trapped at the origin. At the end of inflation the field S starts rolling down the 
potential when H ~ ms- This is the moment of the phase transition, since as S 
approaches zero the effective mass becomes negative. The waterfall field rolls 
down to its minimum <p — v and the effective curvaton mass increases due to its 
coupling with 0. 

The superpotential in Eq. ()28j) is invariant under a discreet symmetry under 
which Qs = Qa — 1 an d Qcj, = 0. This symmetry is chosen to avoid superpotential 
terms of the form 

W = A X S 3 + A 2 <x 3 + A 3 3 + A 4 (x 2 + A 5 S 2 0- (30) 

The first two terms give masses to S and a respectively, and should be very small. 
Remember that H osc ~ max[ms, m a ] should be sufficiently small for low scale infla- 
tion to be possible. The A4 and A5 terms lift the flatness of the a and S direction by 
quartic couplings. Moreover they lead to terms mixing the masses of a and S with 
that of the much heavier field 0; these couplings should be sufficiently small to keep 
the a and S fields light. Finally, the A3 term gives rise to a quartic self interaction 
for 0, which will be important only if A3 > h 2 + k 2 . This coupling does not need to 
be anomalous small. 

Without loss of generality we can take and the couplings k, rj real. The only 
phase left is Arg[ST], which only shows up in the term V 3 8hnSa(f) 2 + h.c. This 
term plays no role in the dynamics, it is zero both before the phase transition when 
= 0, and after when S = 0. We will henceforth neglect this term. Writing 
the potential in terms of real fields (denoted by the same symbol), and with the 
rescaline;s k — > k /y/2, h -> h/y/2, 7] -> 77/V2 and v — > v/y/2, the relevant part of 
the potential becomes: 

V = ^ 2 -v 2 ) 2 + ^ 2 S 2 + ^ 2 a 2 , (31) 

with rj 2 = h 2 + n 2 . The vacuum energy Vq = V(0 = 0) = v 4 7] 2 /2. The effective 
masses for the curvaton and waterfall field are 

m 2 a = (m a )l + (hv) 2 , (32) 
ml = -( V vf + (KSf + (haf, (33) 

The contribution of the ha-term to the 0-mass is sub-dominant at all times. The 
curvaton coupling to standard model (SM) particles is not included explicitly in the 
potential, but is assumed to be suppressed. 

We have chosen for a two-field, hybrid inflation type of potential rather than a 
one-field potential to set off the phase transition, as this offers better perspectives 
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for lowering the bound on H*. Consider a one- field potential, Eqs. (|281 l2Uj) with 
k = 0. Then the field has to be light, as the phase transition starts at H ~ m^, 
when starts rolling away from the origin towards its minimum. Consequently, the 
curvaton mass after the phase transition m a = hv = {h/rf)m^ remains small, unless 
h/rj is large. But a large h means a large decay rate T(cx — > 0), leading to early 
curvaton decay. 

Apart from that, a light 0-field has other problems. First of all, there is no 
reason for it to have a small VEV after inflation instead of sitting at its minimum. 
Secondly, a light field with a negative mass is hard to obtain in SUSY theories. 
Either, the mass can be driven negative by radiative corrections (but this means 
large couplings, and hence, a large decay rate), or the mass term is dominated 
by negative soft contributions (implying a non-minimal Kahler potential in gravity 
mediated SUSY breaking). Thirdly, the light field will fluctuate freely during 
inflation, and Sp^/p^ ~ i?*/0* — which is large and non-Gaussian for small 0*. Its 
contribution to the total density perturbations should be negligible small. There are 
two ways in which can imprint its perturbations on the CMB. If p^ contributes 
significantly to the total energy density, acts itself as a curvaton. This contribution 
to the density perturbations is negligible, less than one percent, if 



The ratio Vt^ is evaluated at the time of the phase transition, when it has its maxi- 
mum value (we omit the subscript "pt" to avoid notational clustering). If the field 
is initially close to the origin 0* ~ 50* ~ H*, this implies < 10~ 7 . Further, can 
play the role of the modulating field in the inhomogeneous reheating scenario [2U| ■ 
Indeed, both the curvaton mass and decay rate are a function of 0. The density 
perturbations produced in this way are Sp/p ~ 50/0d ec , and are negligible only if 
the VEV at the time of curvaton decay is sufficiently large, 0d ec 3> °"*- 

All these constraints and problems are avoided simultaneously by making the 
field heavy during inflation. Of course, the density perturbations of the light field 
5" likewise should be suppressed. But this is done much more easily. It requires 
Ps "C (Tx/S*; as there is no constraint on the VEV during inflation one can simply 
take ^> cr*. The curvaton mass and decay rate are independent of S at the time 
of curvaton decay, and therefore reheating is homogeneous. 



The model parameters have to satisfy several constraints: 

1. If the potential energy density dominates before the phase transition it will 
drive a period of inflation, hybrid inflation for Vo dominance and chaotic in- 
flation for S dominance. A long period of inflation should be avoided as this 
erases all traces of the inflationary period at H*. This is assured if Q^, Q$ < 1 
or S 1 * < Mpi and v 2 r] < H pt M ph 




(34) 



4.2 The bound 
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2. The mass of the waterfall field should be positive and large during inflation, 
so that it is trapped at the origin and does not contribute to the density 
perturbations: kS* > H*, r\v. In addition, the density perturbations produced 
by the flaton field S should be sufficiently small, assured if Qs < ICT^/ct*. 

3. There can only be an improvement with respect to standard scenario if m a = 
hv ^> (mcr)*. This is only possible if the vacuum mass of the waterfall field is 
much larger than the curvaton mass: m^/m^ 1. 

4. There is an hierarchy of scales H^ n < T a < H pt ] the curvaton should decay 
before BBN, but after the phase transition, is the decay rate into SM 
matter, which should exceed curvaton decay into any of the hidden sector 
fields, such as and S. Since m a < after the phase transition, curvaton 
decay into <fi quanta is kinematically inaccessible. However, if the fermionic 
superpartners of <fi are sufficiently light, the curvaton can still decay into them 
and Y a > h 2 m a is needed. 

The parameters f2 pt and fl^, defined in Eqs. (FTTH IMj) . give rise to the equalities 

Pa/^pt = Pcf,/^ = Ptot (35) 

Here = V is the potential energy density stored in the waterfall field before the 
phase transition, p a = {m a o*) 2 is the energy density transferred to the curvaton 
field during the phase transition, and p tot = (if pt M pl ) 2 is the total energy density 
at the time of the phase transition. Further f2 pt < saturated when the energy 
stored in Vq is transferred entirely to the curvaton field. From the above equation 
we can solve 

a* \ V J 

^ f2 pt H pt Mpi ^ 

From the first equality in Eq. (J36|) it seems that direct decay with Q pt ~ — >• 1 
is possible for all H pt , and thus the bound on the Hubble constant during inflation 
can be if* > 10~ 14 GeV, the bound obtained in section 13.2.21 However, this would 
require couplings exceeding unity. Imposing h, rj < 1 together with r\ = y/h 2 + k 2 > 
h 4 , the bound for direct decay with Q pt ~ 1 is 



> ^AJHZMZ > 4 x 1(T 7 GeV, (3S 
n Dt ~i h 



4 If one considers the potential in Eq. as an effective potential, not coming from a super- 
potential of the form Eq. Ij28(l . then there is no relation between the couplings. One can take the 
limit rj < h to try to lower the bound. However, this implies m a > and therefore a large decay 
rate T a = h 2 m,j, which is counter productive. 
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where in the last step we have taken H pt , T a — > fibbn and y/rj/h — > 1. This bound 
can be either found by using Eqs. l2o"|) with 



or from the expression for rj in Eq. (}3Tj) above. As discussed in section 13.2.21 the 
bound for Q pt < 1 is the same as for direct decay with Q pt = 1, as long as the 
curvaton density parameter at the end of the phase transition is large enough and the 
decay rate sufficiently small — see Eqs. f)24l |2~B1 l2Hj) — so that curvaton domination 
occurs before decay, i.e., IV < Q pt H pt . 

The bound on the Hubble constant depends on the coupling values. The lower 
bound above is obtained for h ~ k ~ 1 and A — > Abbn ~ 1CT 11 A grav . Then low 
scale inflation with fi* ~ TeV is possible for (mo-)*, H pt < 1CT 6 GeV, < 10 6 GeV 
and fipt > 1(T 10 . 

It could be argued though, see point (4) above, that the coupling h < A. Taking 
all couplings approximately equal h ~ k ~ A — > Abbn ~ 10~ 9 gives the bound > 
10 -2 GeV. This is still compatible with low scale inflation, but the amount of tuning 
is increased. For example, TeV scale inflation is only possible for (m„) t) H pt < 
5 x 10~ 15 GeV, < 10~ 4 GeV and Q pt > 10~ 6 . Moreover, these parameter choices 
violate Eq. f)27j) : the curvaton scenario is jeopardized by thermal evaporation. 

If the decay rate into SM matter is of gravitational strength curvaton decay 
before nucleosynthesis requires ftpfffpt > (H hhn /M pi y/ 3 a* > 1(T 8 GeV for low scale 
inflation with ~ 10 3 GeV. From Eq. (J24j) we also get an upper bound on the 
time of the phase transition H pt < 10~ 6 /f2 p { 4 GeV. There is only a small window 
of H pt for which the heavy curvaton scenario can work; in addition Q pt > 10~ 3 is 
needed. 

4.3 discussion 

Hubble constants much larger than the theoretical, model independent, lower bound 
if* > 10 -14 GeV are needed in the Higgsed curvaton scenario. The reason is that 
the parameters are not all independent. Tuning both fi^ — > 1 and Qpt/fl^ — > 1 is 
impossible for coupling constants not exceeding unity. 

The model has the same general problems discussed in section 13.31 Late time 
baryogenesis and at least two small masses. More generically there are three small 
masses, but with a distinct hierarchy fi* ^> ^> (m^)*, ms- None of the light 
fields can be a pseudo Nambu-Goldstone boson (PNGB) as the potential terms 
K 2 S 2 <p 2 and h 2 a 2 (p 2 break the shift symmetry of all fields. 




(39) 
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5 The heavy curvaton as axion 



Axions are natural curvaton candidates; the axion mass is protected by a (weakly 
broken) global symmetry, and therefore can be kept light naturally, even during 
inflation when supergravity corrections generically lead to masses of the order of the 
Hubble constant for all scalars. 

5.1 Axion with mass from non-renormalizable operators 

We consider a potential of the form V = Vpq + Vnr with 

V PQ = A 2 (£ 2 - ff. (40) 

The VEV of £ breaks the PQ-symmetry spontaneously at a scale (£) = /. The 
curvaton a is identified with the canonically normalized axion field, defined as £ = 
(//v2) exp(ia/f). If the non-renormalizable potential breaks the PQ symmetry 
explicitly it gives a mass to the axion m a oc f a , with a a model dependent constant. 
Therefore, one way to implement the heavy curvaton is to have an increase in the 
VEV of £ at H pt [HI- This can be obtained by a potential term /i£ 2 2 , with h < 
0. Although negative couplings are possible in SUSY theories, they require very 
elaborate constructions j2H]- We will therefore pursue another option. 

Consider the non-renormalizable potential with operators of dimension D\ = 
2mi + rii and _D 2 = 2Z 2 + 2rri2 + n 2 respectively 

V 1 2m\ 



^nr = A 2 [ M2mi+Wl _ 4 + c.c. + A 2 I + c.c. ) . ( I I ) 



|2Z 2 | S |2m 2 jyi2 
721 

'pi / \ J "V 

This potential gives two contributions to the curvaton mass m 2 = m\ + m\ with 

* - x K£) Dl ~ 2 ^ < 42 » 

/ fD 2 -2h-2j ) 2l 2 \ 

m l = AU J mD J jM 2 . (43) 

During inflation cf) has a small VEV and (m^)* = m\. After inflation there is a phase 
transition in which the VEV of grows to large values and m a ~ m 2 . The bound 
on the inflationary scale can be weakened only if m<i > rri\ or 

•A 2 ' 1 AAV / ~- Dl - D ' 



fj > [i) [tJ • (44) 

For operators of the same dimension D\ = D 2 and couplings Ai, A 2 ~ 0(1) this is 
just the requirement that <fi > f. 
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In the following we will assume that the initial curvaton VEV is of the order of 
its maximal value cr* ~ /; this assumption minimizes the bound on H*. Curvaton 
decay into SM matter is mediated by non-renormalizable operators, so that the 
PQ symmetry is preserved at the normalizable level. The decay rate then is of 
gravitational strength and the bound is given by Eq. 



AMll 3 H*l 3 

Direct decay corresponds to the limit Q pt - > 1 and / — > cr*. The value of H pt cannot 
be lowered/raised arbitrarily in this model. First of all H pt > mi, otherwise the 
bound is increased by a factor (mi/H pt ). Further mi oc H pt , and thus so is the 
decay rate; H pt should be large enough for decay to take place before BBN. Lastly 
> H pt . This together constrains 



max 



rr \ 1/3 
-"bbn 




< fl^ < TTty. (46) 



If Di > 9 the limit H pt — > H^ n is possible and direct decay with fl pt = — > 1 gives 
the theoretical lower bound H* > 10 -14 GeV of section 13.2.21 However, we still have 
to add the explicit form of the potential that is to trigger the phase transition. 
And as for the Higgsed curvaton, we expect that this will give additional constraints 
on the system. We will discuss two cases below focusing on the possibility of low 
scale inflation with H* ~ 10 3 GeV. Note that independently of the potential, the 
lower bound on Q pt from combining Eqs. (|43| l4*E|) is 



-2/3 

n " >ltrl2 lliRSvJ ' (47) 



-12 / \ 



5.1.1 Vfj, quadratic 

To achieve the phase transition we add a quadratic potential for 

V+ = -m^ 2 . (48) 

The phase transition takes place at H pt ~ when the field starts rolling towards 
large VEVs until the potential is lifted by the terms in Vnr. The 0-field is necessarily 
light; its contribution to the density perturbations is negligible small if Eq. is 
satisfied. 

Assume that at large the dominant term in Vnr has 1% > 2, so that the 
potential increases faster than quadratically and there is no runaway behavior, as 
well as n,2 > 2 so that it generates a mass term for the curvaton. This is the 
most advantageous situation for the heavy curvaton, since then all of the potential 
energy is transferred to the curvaton field and fl pt = f^, with Q pt , defined in 
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Eqs. (Unilllll). This gives the relation m 2 f ~ m^cfi = «/ Vt pt H pt M p i, which can be 



solved for £7 



pt- 



P t 



A 2 M 



pi 



h-l { M, 



'pi 



/ 



(49) 



Eqs. fl43l l4T)j) put an upper/lower bound on H pt , and therefore likewise on Q pt . 
Taking A 2 ~ 0(1) we get 

D2-2/2-2 _2_ 



< Apt < 



pi 



2(D 2 -2/ 2 -3) 
2i 2 -3 



/ y V m p1 y PL V / 

- ai 

the case of low scale inflation with if* ~ TeV and / ~ a* = H^/A. 



(50) 



In Table 1 below the values of Q pt are tabulated for various values of D 2 and Z2 f° r 



(D 2 , k) 


(^pt)min 


(^pt)max 


5D 


(6, 2) 








(7, 2) 


6 x 10~ 9 


1 


2 


(8, 2) 


6 x 10~ 3 


1 


18 


(8, 3) 


9 x 10~ 9 


3 x 10~ 10 


2 


(9, 2) 








(9, 3) 


3 x 10~ 6 


1 


6 


(10, 2) 








(10, 3) 


3 x 10~ 2 


1 


40 


(10, 4) 


7 x 10^ 8 


1 x 10" 5 


4 



Table 1 

The heavy curvaton cannot work for (D 2 , l 2 ) = (6,2), (9,2), (10,2) as these 
operators cannot give Q pt > 10~ 12 , see Eq. (|47|). A low value Q pt < 10 -7 , in order 
to suppress the contribution to the density perturbations, is only possible for 
{D 2 , l 2 ) = (7,2), (8,3) and (10,4). The mass = m x < (H pt ) min for all D x > 7, 

in agreement with Eq. ()46|) . 

The dominant contribution to the curvaton mass is provided by a dimension 
D 2 operator with l 2 > 2. In addition there might be operators in Vnr which are 
pure monomials in <p (terms with D 2 = 2l 2 in the language of Eq. ()41Jl ): call the 
dimension of the lowest order such operator D^. If this latter term dominates the 
potential at large 0, then only part of the potential energy stored in is transferred 
to curvaton energy during the phase transition and Q pt < Q^. For = D 2 the 
monomial always dominates. In this case fl pt = Q^V^/Vd^ <C 10~ 12 for all values 
of (D 2 , l 2 ); the heavy curvaton scenario does not work, as follows from Eq. ()47|) . 
The 0- monomial is sub-dominant, and the curvaton scenario works as before, if 

= D 2 + SD with SD > 2. The values of SD for the different combinations of 
(D 2 , l 2 ) are given in Table 1 above. 

The favored operators, those with a low value of 3D, are the same operators that 
give a low value of Q pt . Not surprisingly, these are the operators with the largest 
value of l 2 , maximizing the mass increase (m 2 /mi). 
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5.1.2 V$ quartic 

In this subsection we take the potential responsible for the phase transition to be 
quartic, of the form 

Vj, = -mj0 2 + r] 2 ^. (51) 

If this is all, and H pt ~ as in the previous subsection, there is no gain in 
introducing the quartic term as it can only be harmful and lead to fl pt < fl^. So 
we assume that the potential is embedded in a hybrid inflation type of potential of 
the form Eq. (|31J1 . If this potential originates from the superpotential in Eq. (|28jh 
then the quartic coupling is restricted: r\ > H*/S* > 1CT 15 for ~ TeV, see point 
(2) in section E~2l In this set up it is possible to have m<f, > H pt , which gives more 
freedom in model building. Further, the 0-field is now heavy during inflation and 
does not fluctuate; consequently the upper bound on fi pt of Eq. (J33|) does not apply. 

For large coupling rj > icr 6 ( D2 ~ 4 ), the potential is lifted by the quartic term at 
large and not by the operators in Vnr, and the VEV is = m^/rj. We can solve 
for Q pt = V NR /(H pt M pl ) 2 , giving 

Using the upper/lower bound on the scale of the phase transition, see Eqs. (JUJEBJ), 
we again get an upper/lower bound on Q pt : 

D 2 -2 2(Z 2 -1) 2(D 2 +<2-3) 4l 2 

J_\ l > (H hhB \ 3h t m4> y fj_\ 3-/ 2 ( 3-h 

M P J \M P J \H ptV J pt \M P J \H ptV J 

(53) 

The lower bound is generically small (fi p t) m in *C 1, and can easily be satisfied. 
The whole issue is whether the upper bound can be large enough, Q pt > 10~ 12 for 
TeV scale inflation. The upper bound is maximized in the limit of small H pt (or 
equivalently T a — > H hhn ) and 77. The mass is related to the scale of the phase 
transition by the relation m^/r] = Vt 1 ^ 2 H pt M p \. 

In Table 2 below we have indicated the values of (D2, h) for which Q pt ~ 1 is 
possible. Large fl pt is not possible for Z 2 = 1, independent of the value of D 2 , so we 
left that option out. 
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( n i \ 
[D 2 , h) 


n 1 


a 


(6, 2) 


V 


10 8 - 10 iU 


(7, 2) 


/ 


10° — 10 


(8, 2) 






(8, 3) 


V 


lO" 9 


(9, 2) 






(9, 3) 


si 


10~ 5 


(10, 2) 






(10, 3) 


V 


io~ 8 


(10, 4) 




< io- 8 



Table 2 

Also listed is the boundary value of ft = [va^j H pt rj) / [m^j H pt r]) maiX for which the 
curvaton scenario can work. Here 



1 x 10 1 



4 x 10 2 



j < i ( OptH^, 
20 ( ^pt^ 



1/4 



1/4 



10 



-12 



io- 16 



n 



1/2 



1/2 



1/6 
1/6 



Do = 6, 



D 2 > 6. 



(54) 



For D 2 = 6 the quartic coupling has to be 77 > 10~ 12 for the quartic term to 
dominate over Vnr, whereas for D 2 > 6 the dominant constraint is i] > 10~ 15 from 
the requirement that <ft is heavy during inflation. 

For l 2 = 2, 3 the curvaton scenario can work for ft < ft < 1, with /5 the value 
listed in Table 2. The smaller /3o the larger the parameter space, as there is more 
room for 77, f2 pt , Qa, and/or T a to vary. For l 2 = 2 the exact value of fto depends 
on f2 pt , and there is a /?o-range: the upper bound corresponds to fl pt ~ 1, whereas 
the lower bound is for Q pt = ft. Note that the upper bound on the curvaton density 
parameter after the phase transition is Q pt > fto- 

For l 2 > 4 the upper bound on Q pt is maximized for large 77 and r CT , as follows 
from Eq. (|53|). The value of fto in this case is an upper bound. For example, for 
(D 2 , l 2 ) = (10, 4) one needs ft < 10~ 8 . In this case couplings r\ ~ 1 are possible. 
The lower bound on Q pt of Eq. (}4Tj) can be reached. 

Non-renormalizable terms with = 2l 2 are negligible small, and therefore do 
not alter the above picture, if 7] > A^(0/M p i)^~ 4 ^ 2 or 



Da 



V > \ 



pt 



Da 



Da 



Mi 



pi 



> (10 



-26n D. 



4> 



where in the last step we have taken \a ~ 1, Q pt — > and H ] 



'pt 



10" 



(55) 



(which 



is the lower bound of Eq. ()46j0. The non-renormalizable potential is negligible for 
f] > 10~ 9 , 10~ 13 , 10~ 16 for Da, = 6, 8, 10. Even for Da, = D 2 , only part of parameter 
space is excluded. This is in sharp contrast with the quadratic model discussed in 
the previous subsection, where Da, — D 2 > 2 is needed. 
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5.2 discussion 



The main difference between the Higgsed curvaton and the axionic curvaton is that 
the latter obtains its mass from non-renormalizable operators instead of through 
a direct coupling. The advantage is that the curvaton can be kept light natu- 
rally. Moreover, since the 0cr-coupling is of gravitational strength, there is no strong 
constraint from curvaton decay into <fi or its superpartners. But the use of non- 
renormalizable operators also has it draw backs. There is less freedom in tuning 
couplings. In particular, the curvaton decay rate is of gravitational strength, which 
constrains the model considerably. Further, it is harder to suppress or forbid un- 
wanted operators, as discrete symmetries might be broken by gravity — indeed, the 
non-renormalizable potential violates the PQ symmetry explicitly. This is especially 
a problem for the quadratic potential of section 5.2.1. 

6 Conclusions 

We revisited the bound on the Hubble scale during inflation in the context of the cur- 
vaton scenario. The bound can be weakened in non-standard settings. One specific 
example is the heavy curvaton scenario. In this scenario the curvaton mass increases 
significantly after the end of inflation, (possibly) triggered by a phase transition. We 
reanalyzed the bound in this set up, exhibiting explicitly the dependence on various 
parameters. We take into account the upper bound on the curvaton mass coming 
from both energy conservation and direct decay. A lower bound on the Hubble con- 
stant of if* > 10~ 14 GeV is obtained, in the limit that the mass increase happens 
just before nucleosynthesis. TeV scale inflation requires H pt < 1 GeV. 

We discussed two implementations of the heavy curvaton in detail. The Higgsed 
curvaton receives a mass through a Higgs mechanism, whereas the axionic curvaton 
receives its mass from non-renormalizable operators. Although in theory a Hubble 
constant if* — > 1CT 14 GeV is possible, it turns out that even obtaining TeV scale 
inflation is hard to achieve. The reason is that in these explicit models the various 
parameters cannot all be varied independently. 

Acknowledgments 

The author is supported by the European Union under the RTN contract HPRN- 
CT-2000-00152 Supersymmetry in the Early Universe. 

References 

[1] K. Enqvist and M S. Sloth, Nucl. Phys. B 626, 395 (2002) 
|arXiv:hep-ph/0109214| ; D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002), 



21 



|arXiv:hep-ph/0110002| ; T. Moroi and T. Takaha shi, Phys. Lett . B 522, 215 
(2001) [Erratum-ibid. B 539, 303 (2002)] |arXiv:hep-ph/0110096| . 

[2] K. Dimopoulos and D. H. Lyth, |arXiv:hep-ph/0209180| 

[3] P. Binetruy and M. K. Gaillard, Phys. Rev. D 34 (1986) 3069; K. Kadota and 
E. D. Stewart, JHEP 0307 (2003) 013; L. Randall, M. Soljacic and A. H. Guth, 
Nucl. Phys. B 472 (1996) 377; E. D. Stewart, Phys. Lett. B 391 (1997) 34; 
E. D. Stewart, Phys. Rev. D 56 (1997) 2019; M. Dine and A. Riotto, Phys. Rev. 
Lett. 79 (1997) 2632; M. Bastero-Gil and S. F. King, P hys. Lett. B 423 (1998) 
27; M. Bastero-Gil, V. Di Clemente and S. F. King, |arXiv:hep-ph/02lToTll 
M. Bastero-Gil, V. Di Clemente and S. F. King, |arXiv:hep-ph/02110T2l ~ 

[4] G. Dvali and S. Kachru, |aTXiv:hep-th/0309095; 



[5] N. Kaloper, Phys. Lett. B 583, 1 (2004) |arXiv:hep-ph/0312002 



[6] D. H. Lyth, |arXiv:hep-th/0308lTo| 
[7] M. Postma, |arXiv:astro-ph/0311563| 
[8] M. Giovannini, |arXiv:hep-ph/0310024| 

[9] T. Matsuda, Class. Quant. Grav. 21, Lll (2004) |arX iv:hep-ph/0312058| . 

[10] D. H. Lyth, C Ungarel li and D. Wands, Phys. Rev. D 67, 023503 (2003) 
|arXiv:astro-ph/0208055| . 

[11] D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003) 
|arXiv:astro-ph/0302209| . 

[12] K. Dimopoulos, G. Lazarides, D. Lyth and R. Ruiz de Austri, 
arXiv:hep-ph/0308015 K. Enqvist, A. Mazumdar and M. Postma, Phys. Rev. 



D 67, 121303 (2003) arXiv:astro-ph /0304187|. M. Postma, Phys. Rev. D 67, 
063518 (2003) jarXiv:hep-ph/0212005| . 

[13] A. Kusen ko and M . E. Shaposhnikov, Phys. Lett. B 418, 46 (1998) 
|arXiv:hep-ph/97094"92| . 

[14] P. J. E. Peebles and A. Vilenkin, Phys. Rev. D 59, 063505 (1999) 
|arXiv:astro-ph/9810509| . 

[15] M. Giovannini, Phys. Rev. D 67, 123512 (2003) |arX iv:hep-ph/0301264| . 

[16] M. K. Gaillard, H. Mur ayama and K. A. Olive, Phys. Lett. B 355, 71 (1995), 
|arXiv:hep-ph/95043"07| . 

[17] K. Di mopoulos, D . H. Lyth, A. Notari and A. Riotto, JHEP 0307, 053 (2003) 
|arXiv:hep-ph/03040"50| . 



22 



[18] C. Gordon and A. Lewis, New Astron. Rev. 47, 793 (2003). 

[19] A. Anisimov and M. Dine, Nucl. Phys. B 619, 729 (2001), 
|arXiv:hep-ph/0008058| ; R. Al lahverdi, B. A. Camp bell and J. R. Ellis, 
Nucl. Phys. B 579, 355 (2000), |arXiv:hep-ph/0001122| . 

[20] G. Dvali, A. Gruzinov and M. Zaldarriaga, arXiv:astro-ph/0303591; L. Kofman, 
|arXiv:astro-ph/0303614| 

[21] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 

[22] M. S. Turner, Phys. Rept. 197, 67 (1990). 

[23] E. D. Stewart, Phys. Lett. B 345, 414 (1995) |arXiv:astro-ph/9407040| . 



23 



